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Abstract 
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these flavour- changing partial widths on all the relevant MSSM parameters and also study 
the non-decoupling behaviour of these widths with the SUSY mass parameters. We find that 
these contributions are sizable as compared to the SM ones, and together with the SUSY- 
QCD contributions they can be very efficient as an indirect method in the future search for 
Supersymmetry. 
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1 Introduction 



Rare processes involving Flavour Changing Neutral Currents (FCNC) provide an extremely 
useful tool to investigate new physics beyond the Standard Model (SM). The strong suppres- 
sion of FCNC in the SM is due to the absence of tree-level contributions and the smallness 
of the loop-contributions. The later, being a consequence of the GlM-cancellation mecha- 
nism [1], naturally enhances the sensitivity of Flavour Changing (FC) processes to possible 
non-standard phenomena. Among the various possibilities, the Higgs-mediated FCNC pro- 
cesses involving down-type quarks are particularly interesting since, within the SM, they are 
additionally suppressed by the smallness of the down-type Yukawa couplings. 

The Minimal Supersymmetric Standard Model (MSSM) [2] provides a natural framework 
where such scalar FC interactions could be significant if the soft SUSY-breaking mass terms 
contain some non-diagonal structure in flavour space. In the minimal-flavour- violation scenario 
of the MSSM, squarks are assumed to be aligned with the corresponding quarks. Flavour 
violation in this case originates from the Cabibbo-Kobayashi-Maskawa (CKM) matrix as the 
only source, and proceeds via loop-contributions, in analogy to the SM. Therefore, its size is 
expected to be very small. However, in the more general scenarios that include misalignment 
between the quark and squark sectors, sizeable contributions to FCNC processes are expected 
to occur. It is indeed the case of the radiative corrections from SUSY-loops in the context of 
S-meson physics that come with factors of tan j3, the ratio of the vacuum expectation values 
of the two MSSM Higgs doublets. These tan /3-enhanced SUSY radiative corrections have 
been studied in a number of processes, including — 5° mixing [3, 4] , leptonic B meson 
decays [5-8], and b ^ decays [9-12], and have been found to be large. In order not to 
be in conflict with present experimental data, these in turn imply some restrictions on the 
parameters measuring the size of flavor mixing in the squark sector which, as we have said, is 
mainly produced from quark-squark misalignment [8, 13-15]. 

Other FCNC processes of interest are related to Higgs physics and are also very sensitive 
to supersymmetric quantum effects [16-19]. In particular, the neutral MSSM Higgs-boson 
decays into bs and sb have been proven to be generated quite efficiently from squark-ghiino 
loops if quark-squark misalignment is assumed [16]. These SUSY-QCD loop contributions 
have the peculiarity of being non-vanishing even in the limit of very heavy SUSY particles 
and, in addition, they are enhanced by large tan /3 factors. Such a non-decoupUng behaviour 
of SUSY particles in the Flavor Changing Higgs Decays (FCHD) can be of special interest for 
indirect SUSY searches at future colliders, as the forthcoming LHC and a next e'^e~ linear 
collider, in particular if the SUSY particles turn out to be too heavy to be produced directly. 
The large rates found for the SUSY-QCD contributions to the Higgs partial decay widths into 
bs and sb [16], as well as to the effective FC Higgs couplings to quarks [5-8, 17, 18], are indeed 
quite encouraging. 

In this paper we complete the genuine SUSY quantum effects in FCHD by the computation 
of the SUSY electroweak (SUSY-EW) one-loop contributions to Higgs-boson decays into bs 
and sb. The dominant SUSY-EW radiative corrections to the effective Higgs-boson couplings 
to quarks have been computed recently in the mass-insertion approximation, including both 
misalignment and CKM induce effects [18]. Here we will perform an exact computation of the 
complete SUSY-EW one-loop contributions from squark-chargino and squark-neutralino loops 
to the flavour-nondiagonal decay rates of the three neutral MSSM Higgs bosons and compare 
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them with the SUSY-QCD contributions. Since no approximation is used, our computation is 
vahd for all values of the characteristic parameter measuring the squark-mixing strength and 
for all tan /3 values. Furthermore, we will explore in detail the dependence of these SUSY-EW 
quantum contributions on the MSSM parameters, and we will study their behaviour in the 
large sparticle-mass limit. Both the numerical results and the asymptotic analytical formulas, 
showing the non-decoupling behaviour, can be of particular interest for future Higgs-physics 
studies at next-generation colliders. 

The paper is organized as follows. Section 2 describes the basis for FC interactions in the 
SUSY-EW sector of the MSSM. The computation of the SUSY-EW one-loop contributions 
to the Higgs-fe-s form factors and Higgs decay widths into bs and sb are outlined in section 
3. The numerical analysis of the FCHD rates and a detailed discussion of the dependence on 
the relevant MSSM parameters are included in section 4, and the large sparticle-mass limit is 
studied in section 5. A set of useful and compact formulas valid in this limit is also derived; 
it is listed in the Appendix, together with other details of the computation. 

2 Flavour-changing interactions in the MSSM 

In the MSSM there are two sources of FC phenomena. The first one is common to the Standard 
Model case and is due to mixing in the quark sector. It is produced by the different rotation 
in the d- and w-quark sectors, and its strength is driven by the off-diagonal CKM-matrix 
elements. This mixing produces FC electroweak interaction terms involving charged currents 
and, in particular, supersymmetric electroweak interaction terms of the chargino-quark-squark 
type, which are of interest to this work. The second source of FC phenomena is due to the 
possible misalignment between the rotations that diagonalizc the quark and squark sectors. 
When the squark-mass matrix is expressed in the basis where the squark fields arc parallel to 
the quarks (the super-CKM basis), it is in general non-diagonal in flavour space. This quark- 
squark misalignment produces new FC terms in neutral-current as well as in charged-current 
interactions. In the case of the SUSY-QCD sector, the FC interaction terms involve neutral 
currents of the gluino-quark-squark type, and their effects on FCHD have been studied in 
ref. [16]. Here we focus on the SUSY-EW interaction terms generating FC phenomena, in 
particular on those of the neutralino-quark-squark and chargino-quark-squark type. The 
first one appears exclusively due to quark-squark misalignment, as in the SUSY-QCD case, 
whereas the second one receives contributions from both sources, quark-squark misalignment 
and CKM mixing. 

We assume here that the non-CKM squark mixing is significant only for transitions be- 
tween the third- and second-generation squarks, and that there is only LL mixing, given by a 
similar ansatz as in [20] where it is proportional to the product of the SUSY masses involved. 
This assumption is theoretically well motivated by the radiatively induced flavour off-diagonal 
squark squared-mass entries via RGE from high energies down to the electroweak scale [21]. 
These RGE predict that the flavour changing LL entries scale with the square of the soft-SUSY 
breaking masses, in contrast with the LR (or RL) and the RR entries that scale with one or 
zero powers, respectively. Thus, the hierarchy LL » LR,RL » RR is usually assumed. 
These same estimates also indicate that the LL entry for the mixing between the second and 
third generation squarks is the dominant one due to the larger quark mass factors involved. 
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On the other hand, the LR and RL entries are experimentally more constrained, mainly by 
b s'-f data [14]. With the previous assumption, the squark squared-mass matrices in the 
{cL,CR,iL,iR) and {sL,S]i,bL,bR) basis, respectively, can be written as follows. 
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T|, Qq are the mass, isospin and electric charge of the quark g; mz is the Z boson mass, 
and S\Y = sin6w contains the electroweak mixing angle 6w The relevant MSSM parameters 
in the SUSY-EW sector are, as usual, the soft SUSY-breaking EW gaugino masses. Mi and 
M2, the /i-parameter, the soft SUSY-breaking scalar masses Mq, Mjj, Mfj, and the soft 
SUSY-breaking trilinear parameters, Ag. Owing to the SU{2)l invariance, Mq^^ = Mq ,, 
Mq ft. For simplicity, we have assumed the soft-breaking trilinear matrices to be 



and M, 



Q,t 

diagonal in flavour space. 

In our parametrization (j2.H) . (j2.2j) of flavour mixing in the squark sector, there is only one 
free parameter, Xll, that characterizes the flavour-mixing strength. For the sake of simplicity, 
we assume the same Xll parameter in the i—c and b — s sectors, writing X = Xll from now on 
for a simpler notation. Obviously, the choice A = represents the case of zero flavour mixing. 

In order to diagonalize the two 4x4 squark-mass matrices given above, two 4x4 matrices, 
j^H for the up-tjpe squarks and R^'^^ for the down-type squarks, are needed. Diagonalization 
yields the squark-mass eigenvalues and eigenstates depending on A. This dependence has been 
studied in [16]; typically, two of the eigenvalues are weakly dependent on A, with mass values 
very close to the case with A = 0, and for the other two, one grows with A and the other 
decreases with it. For A = 0, one recovers the usual pairs of physical flavour-diagonal squarks, 
(&i,&2), (si,S2), (^1,^2), and (ci,C2). 

The following paragraph presents the SUSY-EW interaction terms that are responsible for 
flavour-changing neutral Higgs-boson decays. We will write them in the mass-eigenstate basis. 

The chargino-quark-squark interactions are described by the interaction Lagrangian 
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where d can be either a 6 or s quark; the chargino index is j = 1, 2 for the two physical states, 
and the squark index is a = 1,2,3,4, representing the four physical squark states, g denotes 
the SU{2)l gauge coupling, and the coefficients A^iaj ^Raj listed in Appendix A. 
They include the two above commented sources of flavour- changing vertices, quark-squark 
misalignment and CKM mixing. 

The neutralino-quark-squark FC interactions are described by 
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where the neutralino index is a = 1,2,3,4 for the four physical states (the squark index is 
again a = 1, 2, 3, 4), and the expressions for the coefficients Bf^l^, Bjila^ -^L«a; ^^Raa listed 
in Appendix A. Notice that FC effects originate only from quark-squark misalignment in this 
case. 

The chargino-quark-squark and neutralino-quark-squark interactions specified so far in- 
duce flavour-changing Higgs-boson decays like — * bs and sb [H^ = h^, H^, A^) via 
electroweak one-loop contributions involving virtual squarks and charginos or neutralinos. We 
will study those effects in full detail in the forthcoming sections El IH and El 

For completeness, we list also the remaining interaction terms which are not of the flavour- 
changing type, but enter the Higgs-decay matrix elements and are thus relevant for the present 
work. These are the Higgs-quark-quark, Higgs-squark-squark, Higgs-chargino-chargino, and 
Higgs-neutralino-neutralino interactions, reading 
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1,2,3,4), with the coefficients S^^^^, 



S 



P^Rab specified in Appendix A. The couplings gHxUaua 9h^ 



R,q^ Liji Riji ^Labi 

can be found in ref. [16]. 



3 Generating flavour- changing Higgs decays 



We present in this section the computation of the loop-induced flavour- changing neutral Higgs- 
boson decays into second and third generation quarks, — > bs, sb, for = h^, H^, A^. 
We focus on the one-loop contributions following from the SUSY electroweak sector, i.e. on 
squark-chargino/neutralino loops. One-loop contributions from the SUSY-QCD sector have 
been analyzed in [16]. We follow in this paper the notation introduced in [16]. 

For the partial decay widths, the one-loop matrix element for each decay process H qq' , 
in compact form, can be written as follows. 
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Figure 1: Generic one- loop diagrams for squark-chargino/neutralino contributions to the decay 
H ^ bs. X = X~)X° and q = u,d, accordingly, with indices a,P = 1,...,4 for squarks, 
i,j = 1, ...,4 for neutralinos, and i,j = 1,2 for charginos. These diagrams will be refered to 
as ac, be, cc and dc {(In, bjq, cjv, djsi) for charginos (neutralinos), respectively, from left to 
right. 



with the projectors Pl,r = 1/2(1 T Ts)- Fl and are the form factors for each chirality 
projection L and R, respectively. Fl and Fr follow from the explicit calculation of the vertex 
and non-diagonal self-energy diagrams, depicted generically for the bs case in Fig. ^ Similar 
diagrams appear in the sb case. 



For q = b, q' = s, the form factors read explicitly, with the notation F^^j^ = F^^j^[Hx) for 
a given Higgs boson H^, 
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where summation over the various squark and chargino/neutralino indices is to be understood. 
For the 2-point and 3-point integrals, Bq, Cq, Ch, C12, taken from ref. [22], we have introduced 
an abbreviation for the arguments such that for neutralinos, B = BiMfr , m^o , fTT''^-o) , C = 

C{ml, Mjj ^,1711,171"^ ,m|o,m|o) and C = C(m^, M^^, m^, m|o, '"^J^)' similarly for 
charginos but replacing x° 5 by xij and da,i3 by Ua,i3- The expression for the right-handed form 
factor F^^ can be obtained from F^'' by replacing L ^ R in all the terms given in eq. ()3.2|) . 
The definitions of the k and l factors and of Sg^'' can be found in Appendix A. 

Besides the vertex integrals, the form factors contain contributions from the flavour-non- 
diagonal 2-point functions, denoted by S-^ for x = X^,X~- Eq. ()3.2|1 contains the scalar 
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coefficients in the following Lorentz-decomposition of the non-diagonal self-energies, 
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(3.3) 



The corresponding expressions can be found in Appendix A. 

The results of this computation have been obtained in two independent ways, one without 
and one with the support of FeynArts and FormCalc [23], and agreement was found. Thereby, 
the Feynman rules of MSSM vertices with FC effects had to be implemented in FeynArts, 
extending the previous model ffie^. 

The partial widths for the decays bs, sb {H^ = h^, H^, A^) are simple expressions 

in terms of the form factors given above. Assuming that the final states qq' and q'q are 
experimentally not distinguished, the final results for the partial widths are got by adding the 
two individual partial widths, yielding 



T{H^ bs + sb) 



[1- 



+ 



rrib 



3imj,^-mi-mt)iFrFr + F^n'F, 



R 



(x) j-,{x)* 



-QmsmbiFl^'F'j^ 



+ F^^^Ft>] 



(3.4) 



The dependence of the decay rates and branching ratios on the MSSM parameters will be 
discussed in the next section, and the behaviour in the large SUSY mass limit is investigated 
analytically and numerically in section El 



4 Numerical analysis 

Here we numerically estimate the size of the loop-induced FCHD as a function of the MSSM 
parameters and the mixing parameter A. The GUT relations M3 = ag/a s"^ M2 and Mi = 
5/3 s^/c^ M2 are assumed. For the numerical analysis of the FCHD rates, only values of A (in 
the range < A < 1) that lead to physical squark masses above 150 GeV will be considered. 
The present experimental lower mass bounds on the squark masses of the first and second 
squark generation are actually even more stringent that this value [24], but we have chosen 
here this common value of 150 GeV for simplicity and definiteness. Notice that higher values 
of A (or, equivalently of (523)^^ the usual notation of the mass-insertion approximation) 
are disfavored by the present B meson data involving 6 to s transitions, although are not 
definitely excluded [8, 13-15]. Similarly, in view of the present experimental lower bounds 
on the chargino mass [24], we will consider only \n\ values above 90 GeV. A lower limit of 
M2 > 54.8 GeV at 95% CL, when the chargino neutralino and scalar lepton searches are 
combined, is also considered [25]. 

The MSSM parameters needed to determine the partial widths T{Hrc ^ bs + sb), for 
= h^,H^,A^, are the following six quantities, m^, tan/5, /i, M2, Mq, and A, where we 
have chosen, for simplicity, Mq as a common value for the soft SUSY-breaking squark mass 
parameters, Mq = Mq ^ = Mjy ^^^^ = Mf^^^^y and all the various trilinear parameters to be 

^The model file is available on request. 
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Figure 2: r(/i° — > 6s + sb) in GeV as a function of (tan/3, (GeV)). 



universal, At = Ab = Ac = Ag = A. These parameters will be varied over a broad range, 
subject only to our requirements that all the squark masses be heavier than 150 GeV, > 90 
GeV and M2 > 54.8 GeV. In addition, the extra parameter A measuring the FC strength 
will be varied in the range < A < 1, by taking into account the constraints on the squark 
masses. The masses and total decay widths of the Higgs-bosons have been computed using 
the FeynHiggs program version 2.1 beta [26]. 

Figs. El through ini display the numerical results for T{Hx bs + sb) as functions of the 
MSSM parameters, for the specific value A = 0.5. The following default values of the various 
MSSM parameters. 



have been chosen for the figures, to specify those parameters that are not varied in each 
plot. This set of MSSM parameters is in accordance with experimental bounds for the decay 
6 — s> 57 [27], as we checked with the help of the code micrOMEGAs [28], based on leading 
order calculations [11] and some contributions beyond leading order that are important for 
high values of tan/3 [12]. Nevertheless, for illustration of interesting dependences, we explore 
in the following also a wider range of the MSSM parameters for the decay widths. 

In Figs. 121 through m the MSSM parameters have been grouped into pairs (tan/?, niA) and 
(M2, A) in order to visualize the individual dependences of the FCHD widths for each neutral 
Higgs boson. Figs. El and H show T{h^ ^ bs + sb), T{H° ^ bs + sb) and T{A° ^ bs + sb) 
as functions of the pair (tan/3, m^). A common clear behaviour of all three decay widths is 
the increase with tan (5, yielding maximal FC effects at large tan/? values. In the rest of the 
numerical analysis we have chosen tan/5 = 35, as specified in 1)4.11) . Notice that the behaviour 
for the A^ decay is indistinguishable from the case. 

The behaviour with is less uniform. As one can see from Figs.ElandEl the decay widths 
r(if " —>■ bs + sb) and r{A^ —>■ bs + sb) clearly increase with rriA owing to obvious phase space 
effects, i.e. increasing phase space for larger Higgs-boson masses. In contrast, r{h^ — »• bs + sb) 
shows a less obvious dependence on tua- For small values of rriA and large tan/3 we see a 
sharp decrease of the decay width with niA- This is due to the fact that |sina| decreases 



/i = 800 GeV , Mo 
ruA = 400 GeV , M2 



800 GeV, A = 500 GeV, 
300 GeV, tan /5 = 35, 



(4.1) 



8 



50100 50100 
Figure 3: r(if° ^ bs + sb) and r(y4° ^ bs + sb) in GeV as a function of (tan/5, niA (GeV)). 




Figure 4: T{h^ bs + sb) and T{H'^ bs + sb) in GeV as a function of (M2 (GeV), A (GeV)). 

rapidly between = 100 and 130 GeV and this decrease strongly affects the self-energy-like 
diagrams, which are proportional to sin a (see analytical formulas). Since the mass of the 
lightest Higgs boson reaches a constant value and the m^-dependence of sin a is weak for 
large values of rriA, the decay width is very fiat in this range. 

Fig.|3]shows the behaviour of the decay widths with respect to (M2, A). Since the behaviour 
for the decay is indistinguishable from the H'^ case, we do not show the A^ decay dependence 
from now on. Results for the decays are applied to this case directly. Clearly, for the 
heavy Higgs boson, the decay width grows with M2 up to an approximately constant value 
and decreases with the trilinear parameter A. A less obvious behaviour appears for the /i*^ 
case, depending on values of M2, A and the other MSSM parameters. In this case, we found 
a similar behaviour with respect to M2, but in contrast to the previous case we now have an 
increase of the decay width with growing A up to a maximum value and then a decrease with 
this parameter. 

Now we focus on the behaviour with respect to the other MSSM parameters. First, in Fig.^l 
we show the behaviour of the flavour-changing Higgs decays T{Hx bs + sb) {H^ = , H^) 
as functions of the /x parameter for three different values of m^. The shaded regions in these 
figures correspond to the region excluded by LEP bounds on the chargino mass < 90 GeV. 
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We can see that the width for the decay is approximately symmetric under n — /x, 
depending of the ttia values. In contrast, the r{h^ bs + sb) width is more unsymmetric 
with respect to the sign of /i. Note that all decay widths increase with for < 500 GeV, 
then reach a maximum, and finally decrease. Regarding the behaviour at very small /i values, 
we have also found that the widths do not vanish at /i = 0. The origin of this comes entirely 
from contributions driven by electroweak gauge couplings. 



Figure 5: T{Hx bs + sb) {H^ = h^,H^) in GeV as a function of /i (GeV). Lines with 
triangles, boxes, and diamonds correspond to ttia = 200 GeV, tua = 250 GeV, and uia = 400 
GeV respectively. 

Fig. ini shows the behaviour of the h^{\eft panel) and (right panel) decays as functions 
of the common soft SUSY-breaking squark-mass parameter Mq. The region below Mq = 425 
GeV implies too low and hence forbidden values for the squark masses. The decay width 
has a small value for light Mq due to the fact that chargino and neutralino contributions 
have opposite sign there. For higher values of Mq, the neutralino contributions change sign 
and the partial cancellation disappears, therefore, the decay width increases until it reaches a 
maximum and then decreases for heavier squarks. The previously mentioned cancellation for 
small values of Mq is less obvious for the heavy Higgs, the clearly visible effect is the decrease 
due to the growing squark masses. Notice that the decreasing behaviour is slower in this case. 

In the following we study the behaviour of the corresponding FCHD with respect to A. 
The MSSM parameters are again the ones chosen in ()4.1|) . Given this set of parameters, 
the experimental lower bound on the squark masses restricts the allowed range of the mixing 
parameter to < A < 0.93. Figs. [7| and |H1 contain the branching ratios and show the behaviour 
individually for the neutralino contributions (lines with diamonds), the chargino contributions 
(lines with stars), and the total SUSY-EW contributions (lines with boxes) up to the maximal 
allowed value of A. For about A = 0.8, the branching ratio is around 3 x 10~^ for the lightest 
Higgs boson and 1.7 x 10~^ for and A^. We remark that the SM value for this branching 
ratio is several orders of magnitude smaller, yielding Br{HsM — > &s + s6) ~ 4 x 10"^ for 
musM — GeV. A similar result can be extracted from [29]. 

Investigating the contributions from charginos and neutralinos separately, we found that 
the neutralino contribution increases monotonically with A, being exactly zero for A = 0, as 
expected. On the other hand, the contributions from charginos show explicitly the two FC 
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Figure 6: T{Hx ^ bs + sb) as a function of Mq (GeV). Other parameters as in ()4.1|1 . 
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Figure 7: Br{h^, —>■ bs + sb) versus A for the parameter set ()4.H) . Diamonds (stars) denote 
neutrahno (chargino) contributions, boxes the total SUSY-EW contributions. 



sources, CKM and quark-squark misahgnment. In fact, the chargino contribution is different 
from zero for A = 0. The non-zero value at A = is due to CKM mixing which is not present in 
neutrahno (or gluino) loops. The CKM effect and the effect of squark mixing for small A 7^ 
partially cancel each other, leading to a minimum around A = 0.2. This minimum depends 
strongly on the particular input for the SUSY mass parameters. For larger values of A the 
non-CKM flavour-mixing effect dominates and the branching ratio increases, due to a change 
of sign in F^^ after the minimum. For our choice for the SUSY mass parameters, eq. fl4.1|) . 
the total neutrahno contribution to the form factors, Fl and Fr, is negative for all the studied 
values of A, but for the chargino contribution Fr is always positive and Fl changes sign, being 
positive for values of A < 0.2 and negative for the rest of the studied A values. To illustrate 
this more explicitly, we show in FigjHl the branching ratios as function of A for the range of 
small A below 0.4. The constructive interference of both neutrahno and chargino terms lead 
to a multiple enhancement of the individual contributions in the decay rates. 

We notice also that \Fl\ is larger than for both, chargino and neutrahno contributions. 
It is also expected since the dominant FC effect enters through the LL entry of the squark 
mass matrix. Regarding the size of the total chargino/neutralino contributions and which 
diagram is the dominant one, they depend in general on the particular Higgs boson and the 
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Figure 8: Same as in Fig. [3 for < A < 0.4. 



0.4 



Br{h° bs + sb) 



0.002 



. 0015 



0.001 



. 0005 



0.2 0.4 



D.G 0.8 



A 



Figure 9: Br{h^ —>■ bs + sb) as a function of A for = 150 GeV (line with diamonds), 
rriA = 200 GeV (line with stars) and rriA = 250 GeV (line with boxes). 



other MSSM parameters and again on the value of A. In particular, for A = 0.5, we find 
that the total chargino and neutralino contributions are of the same order, with the chargino 
contribution being slightly larger. These comments apply to all three Higgs bosons. Finally, 
for 0.15 < A < 0.4 the total contribution from neutralinos dominates that from charginos, 
also for the three Higgs bosons (see Fig. |HJ. However, the behavior at very small A values 
is opposite to the previous one, and it is in agreement with the results from the effective 
Lagrangian approach which incorporates the mass insertion approximation [18]. 

The results in Fig. [7| and Fig. |S1 are for = 400 GeV. We want to emphasize that the 
decay rates of are much larger for smaller values of (see Fig. ^ and therefore yield 
larger values of the branching ratio Br{h^ —^bs + sb), e.g. Br{h^ — > 6s + s6) ~ 2 x 10"'^ for 
rriA = 150 GeV and A = 0.85. To illustrate the dependence, we show in Fig. El the results 
for this branching ratio as a function of A for tua = 150 , 200 , 250 GeV. Large branching ratios 
are found for small tua and large A values. 

In summary, the FCHD branching ratios that we have found in this section are quite sizable, 
and are in fact some orders of magnitude larger than the corresponding SM rates, but small 
in comparison with the SUSY-QCD contributions computed previously in [16]. There are, 
however, important interference terms, which modify the SUSY-QCD effects remarkably. The 
combined results for the SUSY-EW, SUSY-QCD and the total contributions to the Br^H^ 
bs + sb) Hx = (/i°, H^) are displayed as a function of A in Fig. [TUl Note that the absolute 
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Figure 10: SUSY-EW (Diamonds), SUSY-QCD (stars) and total (boxes) contributions to the 
Br{Hx bs + sb) — {hP, H^) as a function of A for the selected MSSM parameters. 



value for gluino, chargino (above the minimum) and neutralino contributions grow with A 
separately. The SUSY-QCD quantum contributions arc at least one order of magnitude larger 
than the pure SUSY-EW contributions. However, these later contribute with opposite sign, 
providing an important interference effect. 



5 Non-decoupling behaviour of heavy SUSY particles 

In this section we study the non-decoupling behaviour of squarks, charginos and neutralinos 
of the SUSY-EW contributions to FCHD into bs and sb. This non-decoupling behaviour of 
the SUSY particles also happens in the SUSY-QCD contributions [16] and means that the 
FC effects remain non-vanishing even in the most pessimistic scenario of a very heavy SUSY 
spectrum. The motivation to analyze these effects, including both the SUSY-EW and the 
SUSY-QCD contributions, is that they could provide a very efficient way to search for indirect 
SUSY signals at next generation colliders. 

The origin of this non-decoupling behaviour in the SUSY-EW contributions is, as in the 
SUSY-QCD case, the fact that the mass suppression induced by the heavy-particle propagators 
is compensated by the mass parameter factors coming from the interaction vertices, this being 
generic in Higgs-boson physics. It has been previously analyzed at length in the case of 
flavour-preserving MSSM Higgs decays [30-34]. The non-decoupling contributions to effective 
FC Higgs Yukawa couplings have also been studied in the effective- Lagrangian approach for 
the quark sector [17, 18] and the leptonic sector [35]. In this approach, the FC effects are 
encoded in a set of nonholomorphic dimension-four effective operators which appear due to 
SUSY breaking induced from the radiative corrections. It usually considers just the dominant 
contributions that come from large Yukawa couplings and assumes large tan/? values. It is 
simple, but can not easily implement SU{2) x U{1) electroweak symmetry breaking effects 
which are relevant for small values of tan/5. We will use here instead the full diagrammatic 
approach which has the advantage of taking into account all EW loop contributions and is 
valid for all tan/3 values. Since, on the other hand, we are not using the mass-insertion 
approximation, our results are more general, being valid for all values of the FC parameter 
A. We will see that our results converge in the large tan /3 limit and for small A values to the 
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mass insertion approximation results of the effective Lagrangian approach. 

In order to show analytically the non-decoupling behaviour of squarks, charginos and 
neutralinos in the FCHD, we perform a systematic expansion of the form factors involved, 
and hence of the partial widths, in inverse powers of the heavy SUSY masses and look for 
the first term in this expansion. We have considered the simplest hypothesis for the SUSY 
masses where all the soft breaking squark mass parameters, collectively denoted by Mq, the 
fi parameter, the trilinear parameters, collectively denoted by A and the gaugino masses, are 
chosen to be of the same order and much greater than the electroweak scale Mew, 

Ms Mo Mg Ml M2 ^ fi ^ A:^ M Ew , (5.1) 

where Mq = Mq = Mjy = Mj^ and Mg is the gluino mass. The bino and wino soft-breaking 
masses. Mi and M2, are chosen in our numerical evaluations to follow the GUT relations, 

2 2 

Ml = jMg and M2 = rjMg, with 7 = |% and r] = ^. In order to provide more compact 
analytical results that can be easily used for future phenomenological studies, we have also 
analyzed the case of equal SUSY mass parameters, Ms = Mq = Mg = Mi = M2 = fi = A, 
where 7 = 77 = 1. 

In this section, we present the expansions of the SUSY-EW contributions to the form 
factors in inverse powers of Ms and keep just the leading contribution of this expansion by 
considering that all the remaining involved mass scales m//o, uia, mho,mz, mw and niq are of 
order Mew- To this end, we use the results of the expansions of the one- loop functions and of 
the rotation matrices that are given in Appendix B. In the general case of arbitrary 7 and r], 
we find the analytical results for the various contributions coming from char gino- squark and 
neutralino-squark loops that are collected in Appendix C. We next comment these results. 
First, we can see from eqs. ljdljl to (jClTj) that, taking all SUSY mass parameters arbitrarily 
large and of the same order, 0{Ms), the SUSY-EW contributions to the form factors lead to a 
non-zero value. That is, they do not decouple in the large SUSY mass scenario. Regarding the 
ultraviolet behaviour, we have checked that the various contributions shown in formulas ()C.1|) 
to ()C.7|) are separately finite. A discussion on the relative signs and the sizes of chargino and 
neutralino contributions to the form factors, and Fr, is already included in the previous 
section. In general, the approximate analytical results show similar features than the exact 
results. 

In the following we show graphically some of these features. For definiteness, in all these 
plots we choose Mq = M-„ = u = A = Ms, Mi = |4m5, M2 = 4Ms, tan/3 = 35 and, 
in order to simplify the analysis, we fix the Higgs boson masses to the following particular 
values, rriho = 135 GeV, tjiho = 250 GeV and rriA = 250 GeV. The non-decoupling behaviour 
of T{Hx ^ bs + sb), for = (left panel) and = (right panel), as a function of Ms 
is illustrated in Fig. ^2 Here we have fixed A = 0.5. The exact one-loop results (solid lines) 
and the approximate large Ms analytical expansions of eqs. ()C.1|) to ()C.7|) (dashed lines) are 
shown for comparison. The behaviour of T{A^ bs + sb) is practically identical to that of 
r(if° bs + sb) and is not shown for brevity. We can see in these figures that, for large 
values of M5, the exact partial widths tend to a no n- vanishing value, characteristic of the 
non-decoupling behaviour, which is very well described by our asymptotic results. Notice 
that a proper study of the non-decoupling behaviour requires to take the input Higgs boson 
parameters, i.e the masses and the a mixing angle, at their tree level values. However, for 
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Figure 11: Non-decoupling behaviour of T{Hx bs + sb) in GeV with Mq = = A = Mg = 
Ms, Ml = I^Mg and = ^Mg, for H,, = (left panel) and = (right panel) and 
for tan/3 = 35, A = 0.5, niA = = 250 GeV and rriho = 135 GeV. Exact one-loop results 
in solid lines and expansions given in eq. (jdljl in dashed lines are plotted for comparison. 



illustrative purposes in the case of the lightest Higgs boson, we have prefered not to use the less 
phenomenologically appealing tree level value but to take instead an effective input mass value 
of 135 GeV which includes roughly the leading one-loop corrections for tan /3 = 35 , m^o = 250 
GeV and Ms ~ 0{lTeV). Notice also that Ms has been taken up to very large values just 
to illustrate the non-decoupling behaviour, but the convergence of the exact result to the 
asymptotic one is very good already at moderate Msusy values (say > 600 GeV), what makes 
our asymptotic formulas useful for future phenomenological Higgs boson studies. 

The dependence of the total chargino and neutralino contributions on the FC parame- 
ter A are shown in Fig. ^1 Here we have fixed Ms = 1000 GeV. The four lower lines take 
into account the two FC effects, quark mixing from off-diagonal terms in the CKM matrix 
and squark mixing from quark-squark misalignment, and show the comparative size of neu- 
tralino and chargino contributions. The solid (dotted) lower line is the total exact chargino 
(neutralino) contribution. The dashed lower lines are the corresponding approximate results 
given by our asymptotic formulas of Appendix C. First, we see that our asymptotic formulas 
describe extremely well the behavior with the A parameter for the whole studied interval, 
< A < 1, and this is true for both chargino and neutralino contributions. Second, the 
neutralino and chargino contributions behave as expected at very small A values. As has been 
already discussed in section El the chargino contributions are larger than the neutralino ones 
at these small A values, and at A = the first one is non-vanishing whereas the second one 
vanishes. For larger A values the situation is different. Here we find that for 0.2 < A < 0.7 the 
neutralino contribution is larger than the chargino contribution, which has now a minimum at 
about 0.4, and for A > 0.7 again the chargino contribution dominates. As has been mentioned 
in section El the localization of this minimum varies with the choice of MSSM parameters. 
Notice that for moderate and large A values the alternative and more frequently used mass 
insertion approximation fails. It is clear from our plots that for these A values to assume a lin- 
ear behaviour with A in the form factors is not correct and indeed can give a wrong conclusion 
on selecting the dominant contributions. 
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Figure 12: T^H^ — > 6s + sb) as a function of A. Left panel shows the h'^ decay and right panel 
corresponds with the decay. Exact one-loop chargino contribution in sohd hues. The 
largest contribution corresponds to the effect coming only from misalignment; in the other, 
both possible effects are included. The neutralino contribution in dots. The expansions given 
in eq. ()C.ll to lC.7|) in dashed lines are plotted for comparison. The fixed parameters are taken 
as in Fig. ^2 



The two upper lines in Fig. are the chargino contributions for the simplified case of 
VcKM = 1, which means that the only source of FC effects in this case is quark-squark 
misalignment. We see that the contribution vanishes at A = 0, as expected. Again the solid 
line is the exact result and the dashed line is the approximate asymptotic result. As in the 
previous case, we see the excellent agreement of our large Ms results with the exact ones for 
all A values. The interesting feature here is to compare these lines with the previous ones 
and to notice how large is the effect of misalignment as compared to CKM in the chargino 
contributions. Indeed, we see that, for A > 0.2, these two effects enter with different signs 
and give a total chargino contribution that is considerably lower than if we switch off the 
effects from CKM. This reduction of the chargino contribution is what makes finally the two 
contributions, from charginos and neutralinos, to be of comparable size, and therefore to 
neglect the later is not a good approximation. 

Finally, we present the results for the case where all SUSY mass parameters are equal. 
These are very simple formulas that, besides of being useful for future phenomeno logical 
studies, allow us to compare more easily our results with those of the effective Lagrangian 
approach. The separate results for the various contributing diagrams are given in eqs. ()C.9|) 
to ()C.14|) of Appendix C. The total chargino and neutralino contributions to the form factors 
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for equal mass scales are then given by 
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where oew = f;;- The results for F^ are like the previous ones but replacing nib i^t-s, 
nic — >■ nit and taking the complex conjugate. Notice that in eq. ()5.2|) we can see clearly the 
relative importance of the two FC effects, quark-squark misalignment and CKM, which are 
driven respectively by the two functions F{X) and J(A) whose values at the origin are given 
by F{0) = and J(0) = 2. We can also learn from eqs. ()5.2|1 and ()5.3|1 about the relative 
size of the contributions from Yukawa couplings versus those from pure gauge couplings. For 
instance, for A = 0.5, niA = 250 GeV and tan/3 = 35, the contributions from pure gauge 
couplings are about 50% of the ones from Yukawa couphngs and of opposite sign. 

On the other hand, if we keep just the contributions from Yukawa couplings and neglect 
the contributions from pure gauge couplings, only the first term in eq. ()5.2j) . which goes with 
the quark masses, remains. If we now consider the large tan (3 limit and the small A limit of 



the previous result, take Vckm 
get, 

_ _ 



and use the linear approximation where F{X) ~ — |A, 



we 



OiEW 



nib 



4:71 2niy/ cos j3 



tan/^af)* ) A 



(5.4) 



which agrees with the result from the effective-Lagrangian approach [18]. 

Finally, we briefly comment on the interesting limit called decoupling limit where >> 
niEw, which corresponds to a Higgs sector with very heavy Higgs bosons except h^. Notice that 
when X = the common factor, ((T2^''+tan f5 a[^^*) in eqs. ()5.2p and ()5.3|) is (cos a+tan j3 sin a) 
which in this decoupling limit goes to zero. Therefore, the decoupling of the heavy particles 
indeed occur in this case and we recover the SM vanishing value, as expected. 



6 Conclusions 

In this work we have computed the SUSY-EW quantum effects to flavour-changing MSSM 
Higgs-boson decays into bs and sb. We have used the full diagrammatic approach and therefore 
our results are valid for all tan/? values and for all values of the flavour- mixing parameter A. 
We analyzed in full detail the dependence of the FCHD partial widths, with all the relevant 
MSSM parameters and A, and found that they are very sensitive to tan/3, /x, and A. The 
branching ratios grow with both tan/3 and A and reach quite sizable values in comparison 
with the SM ones, in the large tan/3 and A region. For instance, for tan/3 = 35, tua = 400 
GeV and A = 0.8, we found branching ratios of 3 x 10~^ for the and 1.7 x 10^^ for the 
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and H^. These SUSY-EW contributions are subdominant with respect to the SUSY-QCD 
contributions, but they contribute with opposite sign and important interference effects, which 
modify the SUSY-QCD effects remarkably, appear. 

The most interesting feature of these SUSY-EW radiative contributions is their non- 
decouphng behaviour for large values of the SUSY particle masses. We have analyzed this 
behaviour in great detail and found that these SUSY-EW contributions to the FCHD widths, 
as the SUSY-QCD ones, indeed do not vanish for asymptotically large SUSY mass parame- 
ters. We presented a set of analytical asymptotic results that can be of great utility for future 
phenomenological studies. These results are in agreement with the ones of the effective- 
Lagrangian approach in the large tan/3 and small A limit. However, we checked that, for mod- 
erate tan (3 and A values, some of the contributions which are usually neglected in the effective- 
Lagrangian approach can be sizable and, therefore, a realistic estimate of the branching ratios 
should rely better on the full diagrammatic approach. In particular, the contributions driven 
from pure gauge couplings should not be ignored for these moderate values. 

In conclusion, the results presented in this work indicate that a phenomenological detailed 
study of the Higgs boson decays into bs and sb can be very efficient as an indirect method in 
the future search for supersymmetry. 
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Appendix A 

Here are the abbreviations used in the analytical expressions for the form factors F^j^ ()3.2j) . 

^^,x° _ r(6) n(^) p(^) /^'^° - r('') p(^) 

'^Ll ~ ^ Laa^ Rab^ Lab ''LI ~ ^Raa^LjSa 

x,x° _ r)(^) ,^,x° _ 

'^L2 ~ ^Raa^Lab^Rab ''L2 ~ ^ Laa^ RjSa 

x,x° _ riix) Us) x,x° _ pis) 

'^L3 ~ ^Raa^Lab^Lab ''LS ~ ^Laa^LPa 

^^,X° _ rW 

'^LA ~ ^Laa^Rab^Rab 

^^^X° _ r>{b) p.(x) rp(s) 
'^L5 ~ ^Raa^Rab^Lab 

^^.X° _ rW r)W 

'^LQ — ^Laa^Lab^Rab 

^^.x° _ rW r)W T^{s) 

'^Ll ~ ^ Laa^ Lab^ Lab 
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sma \ / cosa 

— cos a I , a^^ — j sin a 
i sin /? J \—i cos 

correspondingly, and N'^-^ = Nai cos 9w + sin 9w, N'^2 — "-^aisin^^F + Na2Cos6w ■ The 
Higgs-squark-squark couplings gH^^qaqa physical basis are given in Appendix A of [16] 

for the up- type squarks and down-type squarks independently. The quantities R^^\ N, 

U and V are the matrices diagonalizing the mass matrices of the up squarks, down squarks, 
neutralinos and charginos, respectively. 
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Appendix B 



In this appendix we give the expressions required to compute the leading contribution to the 
FCHD partial widths in the large SUSY mass limit defined in Sect where Ms ~ Mq ~ Mg ~ 
Ml ~ M2 ~ ~ y4 ^ Mew- For that purpose we first write the values of the squark, chargino 
and neutralino masses and their corresponding rotation matrices and then the formulae for 
the two- and three-point integrals needed. 

• The expressions for the squark masses and rotation matrices, in the limit of large SUSY 
mass parameters and keeping just the leading contribution, are 

M? ^ M2 , M? ^ Ml , M? ^ M2(l - A) , M| ^ Ml{l + A) (B.l) 
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n,24 — 
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-"■31 
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— ''^32 



-fl-lS — -^^14 — -"-21 — -'^'22 — -'^'33 — -'^34 — -'^'41 — -'^'42 — ^ 

^iA-^^t^nf3) (B.2) 

and similar results for R^^^ just replacing b ^ t, s ^ c and tan/? — cot p. 

The expressions for the chargino and neutralino masses and rotation matrices, in the 
limit of large SUSY mass parameters can be found in [36] . To leading order these masses 
are M-- ~ M2, ~ and M^o ~ Mi, M^o ~ M2, M^o = M^o ~ 

In the large SUSY mass limit the two- and three-point one-loop integrals approach 
their corresponding values at zero external momenta, Co,ii,i2('^q, rrij^, rn^r, ml, m|, ml) ~ 
C'o,ii,i2(0, 0, 0; ml, m^, m^) and i?o,i(p^; mf, m^) — -Bo(0; m{, m^), and we write these later 
as, 

Co(0, 0, 0; ml, m\, m^) = i?^) 
Cii(0,0,0;mi,m2,m3) = -^/io(i?i, i?2) 
Ci2(0,0,0;m2,m2,m2) = -l^/,3(i?i, i?^) 

9 

Bo{0; ml, mj) = A - log^ + gi{Ri) 

1 1 77Z^ 

BiiO;ml,ml) = --A + -log^-giiRi)-g2iRi), (B.3) 

where mi, m2, and ma, represent generically the masses of the different particles inside 
the loops, Ri = mi/m2, R2 = mi/m^ and the explicit expressions for the functions 
fi and gi can be found in [30]. For the simplest case, where Ri = R2 = 1, they are 
/i(l, 1) = /io(l, 1) = /i3(l, 1) = 1 and giil) = (?2(1) = 0. 
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Appendix C 



In this appendix we present the expansions of the SUSY-EW contributions to the form factors 
in inverse powers of Af? where Ms ~ Mq ~ Mg ~ Mi ~ M2 ^ fi ^ A ^ Mew ^-nd keeping 
just the leading contribution. The following results are valid for the most general case of 
arbitrary A, 7, 77 and tan/3, where 7 and rj are defined by, Mi = jMg and M2 = rjMg. The 
form factor is. 
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Ml tan^ 6'vi/ 
2M| 18 



'5(A,7) 



(C.5) 
(C.6) 
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9' 



rUb 



Lcff+dj^ IQ^2 2mw cos (3 
2 {A - /itan/3) 



tan^ 9w 
12 



Ml 



+ 3- 



-/i 



F^^^ can be obtained by replacing all (L ^ i?), ^ m^, mc, and taking the 

complex conjugate. The previous result of eq. ()C1|) is valid for all and tan/5 values and 
keep all the involved quark masses, m^, rrif,, rric and m^, different from zero. The different 
functions of A that appear in the above expressions and their behaviour in the — > 1 and 
7 — > 1 limits are given by, 



Ml 



3 



cos/3(Mf -/i2 
1 



■ (Ml cos /3 + /i sin /5)P(A, 7) 



4 cos /3(M| 



(M2 cos /3 + /i sin (3)P{\ vj) 



24 cos /3 [(M| -/i2) 



(M2 sin /3 + /i cos (3) + 



tan^ 
(Mf^ 



(Ml sin /3 + /i cos /3) 



J(A)A 
(C.7) 



F(A) 
J(A) 
5(A,r/) 

G(A,77) 



2 2A A'^ 

-[(A + 1) ln(A + 1) + (A - 1) ln(l - A) - 2A]; Ihn F(A) -- - 

^[(A + 1) ln(A + 1) - (A - 1) ln(l - A)]; lim J(A) ~ 2 + ^ 



/i KT+A, 



yi + A 

7] 

V 



ln(l + A)-ln(l-A)+^i 
ln(l + A)+ln(l-A)+^i 



V 



V 



A^o ^ ' 3 

5(A,1) = F(A), 



G(A,1) = J(A), 



VTTx 

V 

VTTx 

F(A) 



91 
91 



V 



V 



; ^(A,i) 



F(A)A 
2 ' 
J(A)A 



(C. 



In the following we present our results for equal SUSY mass parameters. The contributions 
from charginos are. 



9^ 



rrib 



IQ-k"^ 2mw COS (3 



a 



{x)* 



+ {V^Ly^KM + V;?kmV;1m) {m + AJ(A))] 



a. 



(x) 



\_{y^CKM^CKM ~^ ^CKM^CKm) "^(A) 



+ (V^^mV^cI-m + ^okmK^km) ^(A)] 



(C.9) 
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IQiT^ 2mw cos P 
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^2 + ^1 



8mL sin^ (3 
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CKM ^ CKM"''t 



Tfl ' 

CRM ' CRM C 



:) F{\) + {V^],,y^^,,ml + V^L.^cKM^ 
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+ {Vc 
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, ytb yc 



CKM * CRM 



) (AJ(A) + F(A)(l + tan/?))] 



(C.ll) 



For neutralinos we get, 
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tan^ e 
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